Search for: All records

Abstract Quasiperiodic oscillations (QPOs) have been recently discovered in the short gamma-ray bursts (GRBs) 910711 and 931101B. Their frequencies are consistent with those of the quasiradial and quadrupolar oscillations of binary neutron star (BNS) merger remnants, as obtained in numerical relativity simulations. These simulations reveal quasi-universal relations between the remnant oscillation frequencies and the tidal coupling constant of the binaries. Under the assumption that the observed QPOs are due to these postmerger oscillations, we use the frequency–tide relations in a Bayesian framework to infer the source redshift, as well as the chirp mass and the binary tidal deformability of the BNS progenitors for GRBs 910711 and 931101B. We further use this inference to estimate bounds on the mass–radius relation for neutron stars. By combining the estimates from the two GRBs, we find a 68% credible range $R_{1.4}=12.48_{-0.40}^{+0.41}$km for the radius of a neutron star with massM= 1.4M_⊙, which is one of the tightest bounds to date. 

Abstract Observations of neutron star mergers have the potential to unveil detailed physics of matter and gravity in regimes inaccessible by other experiments. Quantitative comparisons to theory and parameter estimation require nonlinear numerical simulations. However, the detailed physics of energy and momentum transfer between different scales, and the formation and interaction of small scale structures, which can be probed by detectors, are not captured by current simulations. This is where turbulence enters neutron star modelling. This review will outline the theory and current status of turbulence modelling for relativistic neutron star merger simulations. 

Abstract We present the first seconds-long 2D general relativistic neutrino magnetohydrodynamic simulations of accretion-induced collapse (AIC) in rapidly rotating, strongly magnetized white dwarfs (WDs), which might originate as remnants of double-WD mergers. This study examines extreme combinations of magnetic fields and rotation rates, motivated both by the need to address the limitations of 2D axisymmetric simulations and to explore the physics of AIC under rare conditions that, while yet to be observationally confirmed, may be consistent with current theoretical models and account for unusual events. Under these assumptions, our results demonstrate that, if realizable, such systems can generate relativistic jets and neutron-rich outflows with properties consistent with long gamma-ray bursts (LGRBs) accompanied by kilonovae, such as GRB 211211A and GRB 230307A. These findings highlight the potential role of AIC in heavyr-process element production and offer a framework for understanding rare LGRBs associated with kilonova emission. Longer-duration 3D simulations are needed to fully capture magnetic field amplification, resolve instabilities, and determine the fate of the energy retained by the magnetar at the end of the simulations. 

Abstract Understanding the details ofr-process nucleosynthesis in binary neutron star merger (BNSM) ejecta is key to interpreting kilonova observations and identifying the role of BNSMs in the origin of heavy elements. We present a self-consistent, two-dimensional, ray-by-ray radiation-hydrodynamic evolution of BNSM ejecta with an online nuclear network (NN) up to a timescale of days. For the first time, an initial numerical relativity ejecta profile composed of the dynamical component and spiral-wave and disk winds is evolved including detailedr-process reactions and nuclear heating effects. A simple model for the jet energy deposition is also included. Our simulation highlights that the common approach of relating in postprocessing the final nucleosynthesis yields to the initial thermodynamic profile of the ejecta can lead to inaccurate predictions. Moreover, we find that neglecting the details of the radiation-hydrodynamic evolution of the ejecta in nuclear calculations can introduce deviations of up to 1 order of magnitude in the final abundances of several elements, including very light and secondr-process peak elements. The presence of a jet affects element production only in the innermost part of the polar ejecta, and it does not alter the global nucleosynthesis results. Overall, our analysis shows that employing an online NN improves the reliability of nucleosynthesis and kilonova light-curve predictions. 

Context. Multi-messenger observations of binary neutron star mergers can provide information on the neutron star’s equation of state (EOS) above the nuclear saturation density by directly constraining the mass-radius diagram. Aims. We present a Bayesian framework for joint and coherent analyses of multi-messenger binary neutron star signals. As a first application, we analyze the gravitational-wave GW170817 and the kilonova (kN) AT2017gfo data. These results are then combined with the most recent X-ray pulsar analyses of PSR J0030+0451 and PSR J0740+6620 to obtain new EOS constraints. Methods. We extend the bajes infrastructure with a joint likelihood for multiple datasets, support for various semi-analytical kN models, and numerical-relativity (NR)-informed relations for the mass ejecta, as well as a technique to include and marginalize over modeling uncertainties. The analysis of GW170817 used theTEOBResumSeffective-one-body waveform template to model the gravitational-wave signal. The analysis of AT2017gfo used a baseline multicomponent spherically symmetric model for the kN light curves. Various constraints on the mass-radius diagram and neutron star properties were then obtained by resampling over a set of ten million parameterized EOSs, which was built under minimal assumptions (general relativity and causality). Results. We find that a joint and coherent approach improves the inference of the extrinsic parameters (distance) and, among the intrinsic parameters, the mass ratio. The inclusion of NR-informed relations marks a strong improvement over the case in which an agnostic prior is used on the intrinsic parameters. Comparing Bayes factors, we find that the two observations are better explained by the common source hypothesis only by assuming NR-informed relations. These relations break some of the degeneracies in the employed kN models. The EOS inference folding-in PSR J0952-0607 minimum-maximum mass, PSR J0030+0451 and PSR J0740+6620 data constrains, among other quantities, the neutron star radius toR_1.4^TOV= 12.30_{− 0.56}^{+ 0.81}km(R_1.4^TOV= 13.20_{− 0.90}^{+ 0.91}km) and the maximum mass toM_max^TOV= 2.28_{− 0.17}^{+ 0.25}M_⊙(M_max^TOV= 2.32_{− 0.19}^{+ 0.30}M_⊙), where the ST+PDT (PDT-U) analysis of Vinciguerra et al. (2024, ApJ, 961, 62) for PSR J0030+0451 was employed. Hence, the systematics on the PSR J0030+0451 data reduction currently dominate the mass-radius diagram constraints. Conclusions. We conclude that bajes delivers robust analyses in line with other state-of-the-art results in the literature. Strong EOS constraints are provided by pulsars observations, albeit with large systematics in some cases. Current gravitational-wave constraints are compatible with pulsar constraints and can further improve the latter.